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Abstract 

We generalize the classical lifting and recombination scheme for rational and absolute factor¬ 
ization of bivariate polynomials to the case of a critical fiber. We explore different strategies 
for recombinations of the analytic factors, depending on the complexity of the ramification. We 
show that working along a critical fiber leads in some cases to a good theoretical complexity, 
due to the smaller number of analytic factors to recombine. We pay a particular attention to 
the case of polynomials that are non degenerate with respect to their P-adic Newton polytopes. 
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1. Introduction 

Let F £ K[x,y] be a square-free bivariate polynomial of bidegree (d x , d y ) defined over 
a field K. The lifting and recombination scheme for bivariate factorization consists to 
recombine the analytic factors in K[[a;]][j/] of F into the rational factors of F over K. Up 
to our knowledge, this approach led to the best theoretical complexity for factoring dense 
bivariate polynomials, see (20). However, it has only been developed in the case when 
the fiber x = 0 is regular, that is when F(0,y) is separable of degree dy. In this article, 
we generalize it to the case of a critical (non regular) fiber, both for rational and absolute 
factorization issues. A first motivation for this work is that for fields with few elements, 
a regular fiber might not exist. Although working in a well chosen field extension can 
solve this problem (13), this might have a prohibitive cost (3). A second motivation is 
that a critical fiber brings new combinatorial constraints that might speed up the re¬ 
combination process. In particular, the number of absolute analytic factors to recombine 
necessarily decreazes along a critical fiber, due do the presence of ramification. Our main 
result is the existence of a deterministic algorithm that, given the analytic factors of F 
up to a certain precision m, returns the rational factors of F in small polynomial time 
in the total degree. While the regular case requires a precision m = d x + 1 (20), poly¬ 
nomials with highly z-valuated discriminants might need a higher precision for solving 
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recombinations with linear algebra. However, we show that in positive characteristic, 
the precision d x + 1 is always enough to compute the numbers of rational factors and 
that in zero characteristic the precision 2d x is enough to test irreducibility. Moreover, we 
exhibate different combinatorial tricks that allow to solve recombinations with precision 
d x + 1 in many reasonnable situations (Subsection 4.4). The algorithms we develop here 
are not intended to compete in general with actual implementations, but we illustrate 
on some examples that working along a critical fiber improves the complexity at least in 
some particular cases, especially for polynomials that are non degenerate or locally irre¬ 
ducible along the fiber. The strength of our results depends strongly on the complexity 
of analytic factorization, an issue we have not studied here. 

Main result. The prime divisors of F in the rings K[[x]][y] and K[x,y] are respectively 
called analytic and rational factors. The n-truncated analytic factorization of F is the 
data of the residues modulo x n+1 of the irreducible analytic factors of F. Although it 
is a fundamental step of our algorithm, we do not pay attention here to the analytic 
factorization and we introduce the notation C(n) = C(n, F) for the number of arithmetic 
operations over K required for computing the n-truncated analytic factorization of a 
polynomial F eK[i,j/]. When x = 0 is a regular fiber, it’s well known that C ( n ) C 0(nd y ) 
thanks to the multifactor Hensel lifting (14). In general, analytic factorization is more 
tricky and C[n ) is expected to be closely related to the complexity of Puiseux series 
computation. Our main hypothesis on F is the following: 

(H) F is separable with respect to y. 

We can always reduce to hypothesis ( H ) after applying a separable factorization algo¬ 
rithm. For fields with at least d x (2d y + 1) elements, the cost of computing separable 
factorization is 0(d x d y ) by Proposition 8 in (21). This is negligeable when compared to 
all complexity results we obtain here. Hence, hypothesis (FI) might be restrictive for us 
only for fields with few elements. We denote by: 

• p the characteristic of K. 

• s the number of irreducible analytic factors of F in K(x)[y]. We thus have s < d y . 

• q = \y/d\ the integer part of the quotient of the x-adic valuation v of the y-discriminant 
of F with the minimal degree d of the analytic factors. This complexity indicator q 
will be refined in terms of the resultants and the discriminants of the analytic factors 
(see Section 3). 

• u> the universal matrix multiplication exponent (2 < ui < 2.5). 

In all of the sequel, we assume that fast mulplication of polynomials is used. Hence 
two polynomials in K[y] of degrees at most d can be multiplied in softly linear time 0(d). 

Theorem 1. Letm := max(g, d x + l). There exists a deterministic algorithm that, given 
F £ K[x, y] satisfying hypothesis (H), returns its irreducible rational factorization with 
at most 

- 0(md y s ul_1 ) +C(m) arithmetic operations over IK if p = 0 or p > d x (2d y — 1); 

- 0(kmd y s ul ~ 1 ) + 0(k)C(m) arithmetic operations over F p if K = F p fc. 

We have q £ 0(d x d y ) under hypothesis ( H ) and q can reach this order of magnitude 
(Section 3, Example 2.2). If the fiber is regular, then q = 0 in which case our algorithm 
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specializes to that of Lecerf (20), with a complexity 0{d x d y ). For fields with at least 
2 d y — 3 elements, we can always find a fiber over which q < d x + 1, see Remark 3.7. Over 
such a fiber, we get a complexity 0(d x d y ) +C{d x ). The only difference with Lecerf’s 
algorithm is that we need to compute the truncated analytic factorization along a critical 
fiber, a difficulty that is compensed by an expected smaller number s of analytic factors 
to recombine. It’s an open question to know if C(d x ) C 0(d x d y ). An important case is 
that of non degenerate polynomials, for which all edge polynomials of F(x,y — a) have 
simple roots for all a G Pg. In that case, q < d x and s is strictly smaller to the total 
number of lattice points of all edges (see Section 8 for details). In such a case, the analytic 
factorization reduces after some well chosen monomial change of variables to the classical 
Hensel lifting or Newton iteration strategies. A brute force complexity analysis leads in 
that case to C(d x ) C 0(sd x d y ) but we strongly believe that this result is not optimal. 

Example. Suppose given two coprime positive numbers a and b and a field K of charac¬ 
teristic zero or greater or equal to 2a + b. Let 

F(x, y) = ( y a +x b + y a x b )(x a y b + l)((y - 1)“ + x b + x b (y - 1)“) G K[x, y]. 

Then the curve C C P 1 x P 1 defined by F intersects the line x = 0 exactly at the 
points (0,0), (0,1), (0,oo). The Newton diagrams of F at each of the three points are 
constituted of a unique segment with only two lattice points. Hence, F is necessarily non 
degenerate and locally irreducible at each point. In particular, the analytic factors are 
coprime modulo x, and the Hensel lifting strategy leads to C{d x ) = 0(d x d y ). Hence, the 
all rational factorization requires 0(ab) = 0(d x d y ) operations over IK which has to be 
compared to the classical complexity bounds inherent to the choice of a regular fiber, 
namely 0{d x d y s ul ~ l ) with s the number of rational places over a regular fiber (20). Of 
course, the two complexities will be close as soon as s is small. The difference will be 
more remarkable in the absolute case for which a regular fiber imposes s = d y (see here 
after). In that example, we solve recombinations in the absolute case within 0(d x d y ) 
arithmetic operations over K, while working over a regular fiber would lead to 0(d x d y ) 
operations over K (9). 

Locally irreducible polynomials. This example motivates to introduce an important class 
of polynomials for which our approach leads to a good complexity. We say that F is 
locally irreducible along the line x = 0 (resp. absolutely locally irreducible) if the germs 
of curves (C, P) C (Pg, P) defined by F are irreducible over K (resp. over K) at each 
rational place P of the line x = 0, including the place at infinity. For example, F is always 
locally irreducible along a regular fiber. The previous example is also such a polynomial. 

Theorem 2. There exists a deterministic algorithm that, given F G K[i, y] absolutely 
locally irreducible along x = 0, returns its irreducible rational factorization with one 
factorization in K[y] of degree at most d y plus 

- 0{d x dyS u> ~ 1 ) arithmetic operations over K if p = 0 or p > d x (2d y — 1). 

- 0(kd x d y s uJ ~ 1 ) arithmetic operations over F p ifK = and p > d y . 

In the second case, it’s enough to suppose that F is locally irreducible over K. 

Theorem 2 is not a direct application of Theorem 1 since we can have F locally irre¬ 
ducible with q ss d x d y . This is for instance the case when the projective curve defined by 
F is a rational curve with a unique place along x = 0 and smooth outside this place. If 
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F is non degenerate, checking local irreducibility has a negligeable cost (Section 8). In 
general, this is more tricky. However, it has to be noticed that Abbhyankar developed 
in (1) an algorithm for testing local irreducibility of a germ of curve that do not re¬ 
quire blowing-ups or fractional power series (see also (10) for a generalization to positive 
characteristic). The main ingredient is that of approximate roots and the algorithm uses 
almost only resultant computations. Up to our knowledge, no complexity analysis have 
been done yet. 

Counting factors and testing irreducibility. If we rather pay attention to the number of 
factors, it turns out that we need a lower truncation order ( 0{d x ) for fields of positive 
characteristic), leading to a better complexity. We say that K is an absolute field of F 
if it contains the field of definition of the irreducible absolute factors of F, that is, if 
rational and absolute factorizations coincide (as in the previous example). 

Theorem 3. (1) Suppose that p = 0. Then we can test irreducibility of a polynomial 

F £ K[x, j/] satisfying hypothesis ( H ) with one factorization in K[z/] of degree at 
most d y plus 

0(d x d y s U] - 1 )+C{2d x ) 

arithmetic operations over K. 

(2) Suppose that p > 0 or that K is an absolute field of F. We can compute the number 
of rational factors of F £ K[x,z/] satisfying hypothesis ( H ) with one factorization 
in K[y] of degree at most d y plus 

- 0(d x d v s w_1 ) +C(d x ) operations over IK if p = 0 or p > d x (2d v — 1). 

- 0(kd x d y s “ -1 ) + 0{k)C(d x ) operations over F p z/K = F p fc. 

Note that we need not to suppose that IK is an absolute field of F in the case of 
positive characteristic, leading in that case to a much stronger result. Roughly speaking, 
the underlying reason is that the Frobenius gives an efficient test for that an algebraic 
number in K lie in the subfield K (Section 5). 

Absolute factorization. Finally, we apply our results to the problem of absolute factor¬ 
ization, that is factorization over K. Note that rational factorization can be seen as a 
subroutine of absolute factorization. Given F £ K[x, y\ separable with respect to y, we 
represent the absolute factorization of F as a family of pairs 

{(Pi,gi),..., (P t ,qt)} 

where qj £ K[z] is separable, Pj £ K[x,y, z] satisfies deg z Pj < deg qj, the bidegree of 
Pj(x, y, </>) is constant when <f> runs over the roots of qj and 

t 

p{x,y )=n n e ^[*>2/] 

i = l 9j(<^)=0 

is the irreducible factorization of F in K[x, j/]. This representation is not unique. We 
have that t is smaller or equal to the number r of irreducible rational factors, with 
equality if and only if the q ^s are irreducible. In analogy to the rational case, we denote 
by r = X] beg qj the number of irreducible absolute factors of F. We represent the 
absolute analytic factorization of F in K[[x]][y] exactly in the same way, the ring K[x] 
being replaced by K[[x]] (Section 7). We denote by s the number of irreducible analytic 
absolute factors of F , and we introduce C(n) for the complexity of computing the n- 
truncated absolute analytic factorization of F. 
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Theorem 4. Suppose that p = 0 or p > d x (2d y — 1) and let m := max(g, d x + 1). There 
exists a deterministic algorithm that, given F £ K[x, y] satisfying hypothesis ( H), returns 
its absolute factorization with at most 

0{md y s “- 1 + rd x dy) + C(m) C <D(d x d “ +1 ) + C(d x d y ) 

arithmetic operations over IK. We can take m = d x if F is locally absolutely irreducible 
along the fiber x = 0. 

This result has to be compared to (9), Proposition 12, where the authors get complexity 
0(d u+1 + rd x dy) for absolute factorization, where d is the total degree of F. Note that 
in contrast to Theorem 1, we assume here that IK has cardinality greater or equal to 
d x (2d y — 1). If we only pay attention to the number of absolutely irreducible factors, we 
can avoid this hypothesis and we can deal with the only ( d x + l)-truncation order. 

Theorem 5. There exists a deterministic algorithm that, given F £ K[x, ?/] satisfying 
hypothesis (H), returns the number f of irreducible absolute factors of F with at most 

- 0(d x d y s LU ~ 1 ) + C(d x ) operations over K if p = 0 or p > d x (2d y — 1). 

- 0{kd x d y s ul ~ l ) + 0(k)C(d x ) operations over F p if K = F p fc. 

This result has to be compared to (9), Proposition 12, where the authors get complex¬ 
ity 0(d ul+1 ) for computing the number of irreducible absolute factors. As mentionned 
already, the great advantage of our algorithm is that, when working over a regular fiber, 
the number of absolute analytic factors to recombine is always d y , while working over 
critical fibers reduces this number to s < d y . More precisely, if ej and /,; stand respec¬ 
tively for the ramification indices and residue degrees of the s rational places of C over 
K, the difference beewteen s and d y is measured by the formulas 

r r 

s = < d y = '^2e i f i . 

2—1 2=1 

Hence, the more ramified the fiber is, the more we gain during the recombination step. 
Of course, in counterpart, we have to perform analytic factorization along a critical fiber. 

Example. Here is a very simple illustrating example. Suppose for instance that 
F = {y a + V2x b + x b y a ){y a - V2x b + x b y a ) 

for some coprime integers a, b. Then, the curve F = 0 has only one rational place over 
x — 0, with ramification index e = a and residual degree / = 2. Moreover, F is non 
degenerate with respect to its Newton polytope. It follows in particular that m = d x + 1. 
After some monomial change of coordinates, we can apply an absolute Hensel lifting 
strategy which leads to C(d x ) C 0(d x d y ). Since both s and f are constant, it follows 
from Theorem 4 and 5 that we compute the number of absolute factors and the absolute 
factorization of F with respective complexities 0(d x d y ) and (D(d x dy), which have to be 
compared to the complexities 0(d x d y ) inherent to the choice of a regular fiber (9). Of 
course, this is a very special example. In general, it would be really interesting to know 
both in practice and in theory when one approach is better than an other. 
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Main line of the proofs. The approach we propose to solve the problem of recombina¬ 
tions of analytic factors follows closely that of Lecerf (20). Namely, we use logarithmic 
derivatives in order to reduce a multiplicative recombination problem to an additive re¬ 
combination problem. Then, a simple observation shows that we need to test if some 
rational function G/F has all its residues pk s in the subfield K C K(x), a problem that 
can be reduced to a divisibility test by F for zero or big enough characteristic. Note that 
in contrast to (20), we do not make any assumption on the cardinality of the field so 
that we need to take care to the case when the leading coefficient of F is not invertible 
in K[[x]]. For small positive characteristic, we test x-independance of the residues thanks 
to an Fp-linear operator introduced by Niederreiter for univariate factorization (22) and 
extended to the bivariate case by Lecerf (20). Hence, linear algebra over F p appears, 
explaining that our complexity results are expressed only for finite fields when the char¬ 
acteristic is small. When the fiber is regular, residues in IK turn out to be a sufficient 
condition for solving recombinations. Along a critical fiber, this is not the case anymore. 
The basic idea is to introduce extra linear equations that depend on the higher truncated 
analytic factors. To this aim, we introduce the separability order of F which in the rnonic 
case, coincides with the maximal x- valuation of d y F(x, <f>) when </> runs over the roots 
of F. We show that this integer gives an upper bound for the required precision. If we 
know moreover that the residues of G/F lie in the subfield K C K, we show that we can 
improve this upper bound. This is the kind of arguments that allows us to prove (2) in 
Theorem 3. Finally, we extend our results to the absolute case by using a Vandermonde 
matrix that allows to reduce K-linear algebra to K-linear algebra. 

Organization of the paper. In Section 2, we introduce our main notations and we explain 
the recombination problem. In Section 3, we solve the recombination problem along a 
critical fiber. In Section 4, we pay attention to the subproblem of counting the number of 
factors and we give in particular an irreducibility test. We discuss moreover some com¬ 
binatorial approaches for solving recombinations of some so-called reasonnably ramified 
polynomials. In Section 5, we give explicit equations for constant residues, mainly fol¬ 
lowing (20). In Section 6, we develop the algorithms underlying Theorem 1 and 3 and we 
study their complexities. We consider the case of locally irreducible polynomials and we 
prove Theorem 2 in Subsection 6.3. In Section 7, we pay attention to absolute factoriza¬ 
tion and we prove theorems 4 and 5. In Section 8, we consider the case of non degenerate 
polynomials with respect to their P-adic Newton polytopes. We conclude in Section 9. 


2. Factorization, recombinations, residues. 

We explain here the strategy developed by Lecerf in (20) for solving recombinations 
in the regular case, and we show that some problems occur when working along a critical 
fiber. For convenience to the reader, we tried to follow the notations of (20). In all of 
this section, we suppose that F is primitive with respect to y, a situation that can be 
reached with a negligeable cost for our purpose. For convenience, we only pay attention 
to rational factorization, the absolute case being treated separately in Section 7. 
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2.1. The recombination problem 


We normalize F by requiring that its leading coefficient with respect to y has its 
first non zero coefficient equal to 1. The polynomial F thus admits a unique rational 
factorization 

F = F 1 -■■ F r eK[x,y], (1) 

where each Fj £ K[x, y] is irreducible, with leading coefficient with first non zero co¬ 
efficient equal to 1. On another hand, F admits a unique analytic factorization of the 
form 

f’ = «J 1 -J s eK[[ 3 ;]][ ! ,] (2) 

where the Fi £ K[[i]][j/] are irreducible with leading coefficient x n \ rii £ N and u £ K[x], 
u( 0) ^ 0. Hence, each rational factor Fj has a unique normalized factorization 

F j = Cj-Tv* 1 • • • Ff jr , j = (3) 

for some polynomial Cj £ K[x], Cj(0) = 1. The recombination problem consists to compute 
the exponent vectors 

Vj = S N r 

for all j = 1,..., r. Then, the computation of the Fj's follows easily. Since F is squarefree 
by hypothesis, the vectors Vj form a partition of (1,..., 1) of length r. In particular, they 
form up to reordering the reduced echelon basis of the vector subspace they generate 
over any given field F. In positive characteristic, our algorithm will have to solve linear 
equations both over IK and over F p . Hence, in order to unify our notations, we consider 
for a while the recombination problem over F a fixed given subfield of K. Namely, we 
want to compute a basis of the following F-vector space 

S := (v i,... ,v r ) ¥ C F s . 

Hence, solving recombinations essentially reduces to find a system of F-linear equations 
that determine S' C F s . If not specified, all vector spaces we introduce from now are 
defined over F, keeping in mind that F will have to play the role of K or F p . 

Truncated functions. Given Q £ K[[x]][y], we denote by [G] n £ K[x, 2 /] the canonical 
representant of G modulo ( x n ). We call it the n-truncation of Q. We will use also the 
notation 

[Q]n :=G- [G] n and \Q\™ := [G] m - [5]" 

for lower truncation of functions, with convention that [G]™ = 0 for n > m. In other 
words, we put to zero all coefficients of monomials with x-degree < n. 

2.2. Recombination and residues. 

The key point to solve recombinations is to reduce a multiplicative problem to a 
linear algebra problem thanks to the logarithmic derivative operator. Let Fi stands for 
the quotient of F by Fi. Let p = (/zi,..., p s ) €E F s . Applying logarithmic derivative with 
respect to y to (3) and multiplying by F we get the key characterization 

s r 

H £ S ■<=> 3 «i,..., a r £ F | HiFdyFi = ^ ajFjdyFj. (4) 

i=i i=i 
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The reverse implication holds thanks to the separability assumption on F ((20), Lemma 
1). The key idea is to derive from (4) a system of linear equations for S that depends 
only on the ( d x + l)-truncated polynomial 


G ;/ :=Y J ^id y F i 


\d x +l 


[x,y\- 


i= 1 


Let G = Gp and let us denote by pk = pk{p) the residues 

G{x,y k ) ^ iyrr 

Pk ■= n _^ £ K(x), 


d y F(x,y k ) 


k — 1, ... ,dy 


of G/F at the roots y k £ K(x) of F. These residues are well defined thanks to the 
separability assumption on F. We get from (4) that 


p£S ==> pk £ F V A; = 1,..., d y . (5) 

In particular, we have an inclusion of F-vector spaces 

S' C U(L) := eF s | Pk el, k = 1,...,^} 

for any subfield L C K(x). In the regular case, the reverse inclusion holds as soon as 
LcI, thanks to the following proposition ((20), Lemma 2). 


Proposition 2.1. Suppose that F( 0,y) is separable of degree d y . Then S = U(K). 


In characteristic zero or high enough, we get that p £ S if and only if p' k = 0 for all 
k , a condition that can be traduced into a finite number of linear equations over K. In 
small positive characteristic p , we have that p' k = 0 implies that p k £ K(x p ) and we use 
then the Niederreiter operator in order to get some extra F p -linear equations that allow 
to test pk £ K (see Section 3). 

Unfortunately, the equality S = U(K) in Proposition 2.1 no longer holds along a 
critical fiber, as illustrated by the following example. 


Example 2.2. Let F = y 6 — (y — x ) 2 £ Q[x, y]. We see that F(0,y) has a double root 
so that the fiber x = 0 is critical. We compute that F has s = 5 irreducible analytic 
factors over Q and r = 2 rational factors. Two of the analytic factors of F have x-adic 
expansions 

Fi = y — x — x 3 + • • • , F 2 = y — x + x 3 + ■ ■ ■ 

Since d x = 2, it follows that 

[PidyFi} dx+1 = [F 2 d v F 2 \ dx+1 . 

In particular, the vector p = (1,—1,0,0, 0) £ F 5 gives the null polynomial = 0, 
hence the trivial relation p £ V(Q). On another hand, we can check that p ^ S. Hence, 
Proposition 2.1 doesn’t hold in that case. 


This example suggests to quotient V(K) by the vector subspace Z of relations = 0. 
Unfortunately, we could not prove that the isomorphism S — V(K)/Z always hold, 
although this is the kind of approach we will follow in order to compute the number 
of irreducible factors (Subsection 4.1). Moreover, even if such an isomorphism holds, it 







does not allow in general to compute the reduced echelon basis of S thanks to linear 
algebra (see Subsection 4.4). Hence, we rather privilegiate to reduce recombinations to 
linear algebra. To do so, we need to introduce extra equations for S. Not surprisingly, 
these equations will depend now on the analytic factors of F truncated up to some higher 
precision, this precision being closely related to the valuation of the discriminant. 


3. Recombinations along a critical fiber 

In Subsection 3.1, we introduce the notion of separability order of F. This integer will 
measure how much the fiber x = 0 is critical for F and will play the role of an upper 
bound for the truncation order of the analytic factors. In Subsection 3.2, we solve the 
recombination problem along a critical fiber. We keep the same notations and hypothesis 
as in the previous section. In particular, F is primitive with respect to y. 

3.1. The separability order 

To each analytic factor T t of F , we associate the integers 

Ti .— val x Res v (Ji,7i), S% ■ —val x DisCy(.7^), d j .— deg yifTf). 

Here, Res y and Disc y stands for the usual resultants and discriminants with respect to 
y , and val x stands for the x-adic valuation of K[[x]]. We introduce the rational number 


n + Si 



and we denote by N = N(F) the integer: 

N := max { [qi \,..., 

The integer N measures in some sense how critical the fiber x = 0 is for the curve F = 0. 
We call it the separability order of F along the fiber x = 0. In particular, we have N = 0 
if F(0,y) is separable of degree d y (the converse is false, take for instance F = y 2 — x). 
The integer N will play the role of an upper bound for the truncation order that allows 
to solve recombinations. The following lemma summarizes its main properties. We recall 
that the standard x-adic valuation val x of the complete field K((x)) uniquely extends to 
a valuation on its algebraic closure K((x)), that we still denote by val x . 

Lemma 3.1. (1) We have equality 

q\d\ H-1 - q s d s = val x Disc y (F). (6) 

(2) Let (f> be a root of Ti, and denote by rq the x-valuation of the leading coefficient of 
Ti. We have the relation 

qi = val x d y F((j>) + dy . (7) 

CH 

In particular, if the leading coefficient of F is invertible in K[[x]], we have that 
{qi, ■ ■ ■, q s } = { val x d y Fff>), cf roots of T 1 }. 
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Proof. By the multiplicative properties of the discriminant and the resultant, we get that 


Disc y (F) = JjDisCjX-Fj) ReSy(F, F,) 2 

i— 1 

s s 

— Disc y i ) R,6Sy i , P"i ), 

i— 1 

and (1) follows directly by applying vaR to this equality. Let now <f> be a root of F;. We 
thus have 

5 y F(0) = TtdyTi(4). 

On another hand, we have the product formula 

fidyFM) =ReSy(J r i,T i dyJ r i)lc(J r i) 1 ~ dy , 

F ( 0)=0 

where lc(Fj) stands for the leading coefficient of F* and where the left hand side product 
runs over all roots of F,. Combined with the multplicative property of the resultant 

ReSy (J~i , J~idyJ ~,) — ReSy (J ~\, F*>) ReSy ( F <, S) y J ~\) 

and with its relation to the discriminant 

ReSy ( J~ i . OyJ~ j ) lc(Fj) DiSCy(F), 

we get the formula 

TidyT^) = lc (Ti) 2 ~ dy ReSy(Fi, F) Disc y (FR. 

f(W)=o 

Since vaR is invariant under the K((x))-automorphisms of the algebraic closure of K((x)), 
point (2) follows by applying vaR to the previous equality and by dividing by the degree 
di of F;. □ 

Remark 3.2. Lemma 3.1 implies in particular that 

^ vaR DisCy(F) ^ d x (2d y - 1) 
d ' d 

where d := min{<R, i = 1,..., s} stands for the minimal degree of the F’s. In particular, 
N £ 0(d x dy). The following generalization of Example 2.2, suggested to us by Eduardo 
Casas-Alvero, shows that N may reach this order of magnitude. 

Example 3.3. Let F(x, y) := (y — x m ) 2 + y n £ Q[x, y], with n > 3 odd. Then (0, 0) is 
the unique point of the curve F = 0 that is ramified over x = 0. We can show that F 
admits a unique irreducible analytic factor T\ vanishing at (0, 0), with degree d\ = 2. It 
follows that 

R = vaR DisCy(F) and r\ = 0, 

while Si = ri = 0 for alH > 1. We compute here that vaR Disc y (F) = mn. It follows that 


n + | 

mn 


dx dy 

L di J 

l 2 J 


l 4 J 


which is of the order of magnitude of d x d y . 
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3.2. Solving recombinations along a critical fiber. 


We can derive from (4) an other obvious source of equations for S. Namely, let us 
introduce for n £ N the F-vector subspace 

S 

W n := [p £ F s I = o}, 

i=l 

with convention W” = F s when n < <4, + 1. For a question of degree, (4) implies that 
we have the inclusions 

ScW n VneN. 

Our next result ensures that the separability order gives an a priori upper bound for n 
for which W n provides enough extra equations to solve the recombination problem. 

Theorem 3.4. We have S = T(K) D W n for all n > N. 

In particular, if TV < d x , then the recombinations are solved by the same system of 
linear equations as in the regular case: 

Corollary 3.5. Suppose that N <d x . Then S = V(K). 

Remark 3.6. This corollary implies in particular that all polynomials that are non 
degenerate with respect to their Newton polytope satisfy V(K) = S (see Section 8). 

Remark 3.7. Suppose that F is separable with respect to y. For a £ Pjjj, let us denote 
by N a the separability order of F over the fiber x = a. From inequalities, 

^2 N a < ^2 vala,_ a Disc y (F) = deg^, Disc,, (F) < d x (2d y - 1), 

we deduce that there always exist a fiber for which N a < d x as soon as K has cardinality 
> 2 d v - 3. 

In order to prove Theorem 3.4, we need to prove two preliminary lemmas. The first 
lemma is a key lemma that will be used many times in the paper. 

Lemma 3.8. Let KcLc K((x)) be a field and let G £ K[x,y], with deg y G < d y . The 
residues pk of G/F all lie in L if and only if G is h-linear combination of the Ejd y Ej’s, 
where the Ej’s stand for the irreducible factors of F over L. 

Proof. One direction is clear: if G = 'YhotjEjdyEj, then pk = ay- £ L where j is de¬ 
termined by condition Ej(x 7 yk) = 0. Suppose now that pk £ L. Thanks to the degree 
assumption on G and the separability assumption on F , we have the partial fraction 
decomposition 

d y 

\ ' pk 

r~tiy~y k ' 

Let r £ r := Aut ( L(x)/L(x)) acts on this equality. By assumption, r leaves both G/F 
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Pk 


and pk fixed. Hence, we get 

dy 

E 


Pk 


*;=i 


y-T{yk) 


= E^- = E 


Pk r 


fc=1 


2/ — 2/fe 


fe=i 


2/ - T(yfc) ’ 


the second equality using that r permutes the roots of F. Here, the notation k T stands 
for the unique index such that (fk T = r(j/fc). The partial fraction decomposition being 
unique, previous equality implies that 

Pk = Pk T V r e T. 

Since T acts transitively on the set of roots of each L-irreducible factor Ej , it follows that 
Pk = Pk’ whenever Ej ( x , yk) = Ej (x, yk' ) = 0. Hence, there exist constants Qi,...,a<6L 
such that 

t 


G 

F 


S>>( E 


1 


y-yk 


3 =1 k\Ej (y k )=0 

The result follows from multiplication by F. □ 


= E° 

3=1 


®yEj 
° Ej ' 


The next lemma computes the valuations of the roots of F. 


Lemma 3.9. Let F £ K[[x]][y] be an irreducible polynomial of degree d in y. Let a and b 
stand respectively for the valuation of the leading coefficient and the constant coefficient 
of F seen as a polynomial in y. Let <j> £ K((x)) be a root of F. Then val c (<j>) = (b — a)/d 
and either a or b is equal to 0. 


Proof. Since F is irreducible, at least one of its coefficient has valuation 0. Hence, if 
both a and b are non zero, then its Newton polytope would contain at least two distinct 
compact edges (Section 8). This is impossible since F is irreducible. Let N stands for the 
norm of the field extension of K((*)) defined by F. Then N(</>) is equal to the quotient 
of the constant coefficient of F by its leading coefficient. Hence val^ N((/>) = b — a and we 
conclude thanks to the relation val x </> = val x N(</>)/ deg(</>). n 

Proof of Theorem 3.f. We already saw that S C V(K) fl W n and we need to prove the 
reverse inclusion when n > N. Let y £ V(K) fl W n . Thanks to the previous lemma, and 
by definition of W n , we deduce that there exists some constants a.j £ IK such that 

S i 

E W [M -^] m = E ^EjdyEj , (8) 

i=1 3=1 

where m = max(d 2 ; +1, n) and where the Ej's stand for the irreducible factors of F over 
K. Let £ K((x)) be a root of Fi and let j be the unique index such that Ej{cj)) = 0. 
Using the relations 

FidyFi{(j )) = EjdyEj((f>) = d y F ((/>), 
we get by evaluating (8) at <f> an equality 

(Pi - aj)8 V F (</») = x m R(cf) (9) 

for some R £ K[[x]][y]. We need a lower bound on the valuation of R((f>). We remark that 
the coefficient of y dy ~ 1 in d y F is equal to d y lc y (F). Since the leading coefficient of F is a 
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polynomial in x of degree at most d Xl equation (9) implies that R has y-degree < d y — 2. 
Hence, ultrametric inequality combined with Lemma 3.9 gives 

(d — 2 ^ 77/ ■ 

val x R(<j)) > min{val x 4> l , i = 0, ..., d y — 2} > ---, 

di 

(recall that x Ui stands for the leading coefficient of F,). Suppose that /q ^ ay. Hence, 
(9) gives 

val x dyF(4>) > m — ^ v . 

di 

By Lemma 3.1, this is equivalent to that m < qi, contradicting our hypothesis m = 
max(d x + l,n) > N. It follows that yt = ctj. Combined with (8), we get that 


G := 


S 

^ ^ ( \llG~ 


d x +1 


s 

^ [dyJ~i . 

i-1 


In particular, the residues of G/F all lie in Fcl, and it follows from Lemma 3.8 that 


r r 

G = 'y ' HiTjdyJ-i = 'y ~ CjFjdyFj. 

i=s i=1 


for some Cj £ K. Since Jq is coprime to Fid y Ti by hypothesis, this relation forces equality 
yi = Cj when Ti divides Fj. It follows that y £ S. □ 


4. Counting the number of irreducible factors 

We show here how to bound the number of factors with the d x + 1-truncation order 
and we deduce a deterministic irreducibility test that requires the only 2d x -truncation 
order. We still suppose that F is primitive with respect to y. 

4.1. An upper bound for the number of factors 

Example 2.2 suggests to introduce the vector subspace Z of vectors y whose associated 
trunacted polynomial Gy, is null, that is 

S 

z := {y e F s I = °}- 

i=l 

We have the following result. 

Proposition 4.1. We have H(F) = S ® Z. 

Proof. We already saw that S C H(F), while the inclusion Z C V(F) trivially holds. 
Hence, we get an inclusion S + Z C H(F). Let us show the reverse inclusion. If y £ V (F), 
it follows from Lemma 3.8 that Gy is F-linear combinations of the irreducible factors of 
F over K. It follows that 

S S 

Ytk[F t dyF l } d ^ +1 = ]Ta i [.F l S y j;] d * +1 , 

2=1 2=1 


13 





for some a = (op,..., a s ) £ S. In particular, /i — a £ Z. Equality V(F) = S + Z follows. 
Finally, if p £ S D Z we get that 

S S 

[liTidyTi = ^ M* [TidyJ r i) dx + 1 = 0, 

i=l t=l 

so that fj, = 0 by linear independance of the FjOyFfs. It follows that V (F) = S ® W. □ 

Remark 4.2. We have Z = 0 as soon as the separability order satisfies N < d x + 1. 
Namely, we have in that case S = V(K) by Theorem 3.4 and we conclude thanks to the 
inclusion S ® Z = V(K) C V(K). 

For fields of positive characteristic, we can take F as the prime field of K, in which case 
the Niederreister operator leads to an explicit system equations for V(F) (see Section 5). 
Hence, Proposition 4.1 allows to compute the number of irreducible factors 

r = dimp V (F) — dimp Z 

with linear algebra from the (d x + l)-truncated analytic factors only. For fields of char¬ 
acteristic zero, testing whether the residues lie in K is a much harder task. In that case, 
we only get equations for V(K), so that Proposition 4.1 a priori allows only to compute 
the upper bound 

r < dimp V (K) — dimp Z. 

This problem motivates to explore in more details the relations beetween V (F) and V (K). 

4-2. On the relations beetween V(F) and V(K) 

In regards to the Proposition 4.1, we may ask whether equality P(F) = V(K) holds. 
We could not prove nor disprove this equality. However, we give here some conditions 
under which it holds. Let us first note the following lemma. 

Lemma 4.3. If all the absolute factors of F are defined over K, then H(K) = V’(K). 

Proof. By Lemma 3.8, if ji £ V(K), then ^ cijEjd y Ej for some ctj £ K, and 

where the Ej 's stand for the irreducible absolute factors of F. By assumption, we have 
that Ej £ K[x]. Applying r £ AutjK(K) to the previous equality, and using that G M has 
coefficients in K, we get that 

T ( a j)Ej d y Ej = yy a jEjd y Ej, 

3 3 

which implies that r(aj) = otj by K-linear independance of the Ejd y Ej’ s. This being 
true for all r, it follows that aij £ K. Hence fi £ V(1K) by Lemma 3.8. □ 

To each rational factor Fj of F, we associate the integer 

Mj := min|L ( ?iJ> Fi divides Fj in K[[x]][y]|. 

Roughly speaking, Mj measures the minimal contact order of the curve Fj = 0 with the 
complementary curve Fj = 0 along the line x = 0. We denote by 

M := max{Mj , j = 1,..., r}. 

Note the obvious relation M < N with the separability order. 
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Proposition 4.4. We have V(K) D W n = P(F) D W n for all n > M. 

Proof. The proof is similar to that of Theorem 3.4. Let n > M and denote by to := 
max{n, d x + 1}. By Lemma 3.8 and by definition of W n we have y £ V (K) D W n if and 
only if 

S t 

Y HilfidyFiT = J2 a kE k d y E k , (10) 

2=1 k= 1 

where the E k £ K[x, y\ stand for the absolutely irreducible factors of F. Let us fix Fj 
a rational factor of F. By assumption, there exists J~i a divisor of Fj such that < n. 
Let E k be a divisor of Fj. Then E k shares at least one root f> £ K(x) with T r . Hence, by 
evaluating (10) at (f> we get that 

{fH - u k )d y F(4) = x m R{(j)) 

for some R £ K[[x]][y]. Taking a;-valuations, and reasonning as in the proof of Theorem 
3.4, we get that pi ^ a k implies <?i > to, a contradiction. Hence a k = Pi £ F for all 
irreducible factors E k of Fj. Repeating this reasonning for all factors Fj of F, we deduce 
by regrouping the factors E k by conjugacy classes that we have 

S S 

= J 2 cjFjdyFj, 
i= 1 j=1 

for some Cj’s in F. For a degree reason, this is equivalent to that 

s s s 

Y Hi[Fid v Fi] d * +1 - Y c i P i d v F i and E wte]£ +1 = 0 

2—1 j= 1 2=1 

The first equation is equivalent to that p £ P(F) by Lemma 3.8, while second equation 
is equivalent to that /i £ W n by definition. □ 

Last Proposition says in particular that if each irreducible rational factor of F has at 
least one branch with Qi < d Xl then we have P(K) = P(F). This is the case for instance 
in Example 2.2, Section 2. Here is a trivial example that illustrates that the converse 
doesn’t hold. 

Example 4.5. Let F{x,y) = {{y — x ) 2 + y w ){y — x) £ Q[x, y]. Then F has exactly 3 
anaytic factors over Q which satisfy 

T\ = (y — x ) 2 mod x 4 , F 2 = y — x mod x 4 , F 3 = y s + 1 + • • • . 

We find here that q\ = q 2 = 10. In particular, the rational factor y — xofF has a unique 
branch and this branch satisfies q 2 > d x + 1 = 4. On another hand, we have that 

F x = ff mod x d " +1 => (1, -2, 0) £ Z 

and we can show that this is the only possible relation. Hence AivaZ = 1. Since clearly 
dim S' = 2, it follows from Proposition 4.1 that dim V(Q) = 3. Since s = 3 is the 
dimension of the ambient space, it follows that P(Q) = P(Q) = Q 3 . Observe that we 
could not use directly Lemma 4.3 to show this equality since F has two absolute factors 
y — x + iy 5 and y — x + iy 5 that are not defined over Q. 
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4-3. Number of factors. Irreducibility test. 


Proposition 4.4 leads to a formula for r that depends only on the M-truncated factors: 

Corollary 4.6. The number of rational factors is equal to 

r = dim V(K) n W M+1 - dim Z (~l W M+1 , 

hence can be computed with the only truncated precision ma ~x.(d x + 1, M + 1). 

Proof. We know from Proposition 4.1 that V(F) = S © Z. Intersecting with W n , and 
using that S C W n , we get that 

s®(z n w n ) = {s® z) nw n = P(F) nr = v(K) n w n . 

for all n > M, the last equality thanks to Proposition 4.4. The corollary follows by 
counting dimensions. □ 

Of course, we can not a priori compute M without knowing the rational factorization 
so that Corollary 4.6 seems to be useless from a computational point of view. However, 
it leads to an irreducibility test over IK with the only 2d x -truncated precision. 

Corollary 4.7. The polynomial F is irreducible over K if and only if 
dim V{K) O W 2dx - dim Z 0 W 2dx = 1. 

Proof. Suppose that dim V (K) D W 2dx — dim Z 0 W 2dx = 1. Since the inclusion W 2dx C 
W n holds for all n > 2 d x we deduce from Corollary 4.6 that dim S = 1. Suppose now that 
F is irreducible over K. Then it’s enough to show that M < 2 d x by the same argument. 
Suppose on the contrary that M > 2 d x . Since F is irreducible, we have by definition of 
M that qi > 2 d x for all i. It follows from Lemma 3.1 that 

S S 

^^la:(DisCy F'j — ^ ^ Qidi 2 d x ^ ' di — 2d x d y , 

which is impossible for a degree reason. □ 

4-4- ^ combinatorial approach for solving recombinations 

We show here that under some reasonnable conditions, we can compute the factor¬ 
ization of F just by knowing V (F) and Z , hence from the {d x + l)-truncated analytic 
factors only. Let us introduce the subset 

I := ji G {1,..., s}, p G Z => pi = 0 j. 

and let 

L := j//S F s , 

We denote by n : F s —> L the natural projection on L. 

Definition 4.8. We say that F is reasonnably ramified over x = 0 if 

dim n(V (F)) = dim V (F) — dim Z. 
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In other words F is reasonnably ramified if and only if for all j = 1 ,... ,r, there is an 
analytic factor F, of Fj such that p € Z implies pi = 0. In particular, if M < d x + 1, 
then F is reasonnably ramified thanks to the proof of Proposition 4.4. Note that for fields 
with positive characteristic, we can test if F is reasonnably ramified since we can then 
compute 14(F), Z and L (see Subsection 5.2). In characteristic zero, this will be the case 
if we know moreover that V(1K) = 14(F). 


Proposition 4.9. Suppose that IK has characteristic zero or strictly greater than d y . 
Suppose that F is reasonnably ramified over x = 0. Suppose that V (F) = Z © T for some 
vector subspace T whose reduced echelon basis (w±,... ,w r ) form a partition of (1,..., 1). 
Then , 


Fj =prim[lc(P)PJj r “ 


d x -\~ 1 


3 = 


and all rational factors Fj of F can be computed within 0(d x d y ) field operations. 


Here, prim stands for the primitive part with respect to y. In order to prove Proposition 
4.9, we first need a key lemma. We denote by R := K[[a;]]/(a;' i:r+1 ) and by fi the class of 
Fi in R. Since F is not divisible by x, we have that /, £ R*[y], where R* stands for the 
multiplicative group of non zero divisors. In particular, it makes sense to compute // in 
the total ring of fractions of R[y] for any integer k £ Z. 


Lemma 4.10. Suppose that K has characteristic zero or strictly greater than d v , and let 
p £ {0,1} S . Then p £ Z if and only if J{} = l ft' = 1 £ R- If IK has characteristic zero, 
the same conclusion holds with the weaker hypothesis p £ Z s . 


Proof. We have equaliity 

..ins c ( II ft ') 

/ . TifiOyfi ~ f j-j y/Xi 

in the total ring of fractions of R[y). Hence p £ Z if and only if ( fi H ) ' = 0, which is 
equivalent to that J~[ /fi 4 £ R*{y p ), where p stands for the characteristic of K. If p > d y , 
and since p t £ {0,1}, we necessarily have /fi 4 £ R* for a degree reason. If p = 0 the 
same holds obviously. In particular, we must have 

n/r=n^M ; =n ic ^) w= 

for some k £ Z. But we know that val x (fi) = 0 for all i, hence we must have k = 0. □ 


Proof of Proposition f.9. we have by assumption that I4(K) = S ® Z where tt(Z) = 0 
and where the reduced echelon basis (vi,..., v s ) of S is such that ... n(v s ) are non 

zero vectors of {0,1} S in reduced echelon form. Hence the same property has to hold 
for the basis (uq,..., w s ) of T and up to reordering the -Wj's, we must have equality 
7T (wj) = 7T (vj), forcing relations Wj — Vj £ Z. Then the proof of Proposition 4.9 then 
follows from Lemma 4.10 combined with relations (2) and (3). Since Wj has entries in 
{0,1}, the complexity for computing Fj belongs to (D(d x deg y Fj) using fast multiplication 
in R[y] (see the proof of Proposition 6.1 in Section 6 for details concerning complexity 
issues). The last statement then follows by adding this cost over all j. □. 
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Example 4.11. Let us return to example 2.2 of Subsection 2.2. We have that 

V(K) = <(1,0,0,0,1), (0,1,0,0,1), (0,0,1,1,-1)) and Z = <(1,-1,0,0,0)). 

Hence 7 t(E(K)) = ((1, 0,0, 0,1), (0,1, 0,0, 0)) has dimension dim U(K) — dim Z = 2(= 
dim S) and F is reasonnably ramified. Let T be such that V (F) = Z ® T and such 
that the reduced echelon basis {w\ . w 2 ) of T is a partition of (1,..., 1). The constraints 
Wi € {0, l} 5 D V(K), Wi ^ 0 lead to the two possible solutions 

wy = (1,0, 0,0,1) or w\ = (1, 0,1,1,0). 

Using the relation w\ + w 2 = (1,1,1,1,1), these solutions correspond respectively to the 
factorizations 

(Fi,F 2 ) = ([F 1 F 5 ] d ^+\ [F 2 F 3 F 4 \ dx+1 ) or (. F U F 2 ) = ([F 1 F 3 F 4 ] dx+1 , [F 2 F 5 } dx+1 ) 

But we know here that [F\\ dx+1 = [F 2 ] dx+1 since (1, —1,0,0, 0) £ Z. Hence both solutions 
determine the irreducible factorization of F, as predicted by Proposition 4.9. 

Example 4.12. Let us return to example 4.5. We have 

V(K) = <(1,0,0), (0,1,0), (0,0,1)) and Z = <(1, -2,0)). 

We have that 7 t(V(K)) = ((0,0,1)) has dimension 1 < 2 = dimU(K) — dimZ. So F is 
not reasonnably ramified. The family of all possible complementary subspaces T of Z in 
V(K) whose basis form a partition of (1,1,1) is 

T= <(1,0,0), (0,1,1)), T= ((0,1,0), (1,0,1)) or T = ((0, 0,1), (1,1, 0)). 

Contrarly to the previous example, only the second solution leads to the good factoriza¬ 
tion of F. 


Unfortunately, for more complicated examples, looking for a complementary vector 
space T of Z in V(K) whose reduced echelon basis form a partition of (1,..., 1) might 
not be an easy task, even though we know such a T exists. 

An alternative approach in the zero characteristic case is to use linear algebra over Z. 
Namely, we can suppose in that case that F = Q, so that V (F) D Z s is a free Z-module or 
rank dimU(F). Recall that the Hermite normal form of a matrix with integer entries is 
such that the leading entry (first non zero entry) of a nonzero row is positive and strictly 
to the right of the leading entry of the row above it. Moreover, all entries in a column 
above a leading entry are nonnegative and strictly smaller than the leading entry. This 
forces also all entries in a column below a leading entry to be zero. Such a form exists 
and is unique, and it conserves the row space (26). We have: 


Proposition 4.13. Suppose that IK has characteristic zero and that F is reasonnably 
ramified. Then we can order the set {1,..., s} such that I = {1,... ,£} for some i > s. 
Let (w 1 ,..., w r ) be the first r vectors of the Hermite normal form of a basis of the free 
Z-module V(F)nZ s . Then, 


Fj = prim 


ic(F)n^r i 

i=1 


d x +1 


j = l,...,r. 


If moreover the Wj’s have positive entries, then we can compute the Fj’s within 0(d x d y ) 
field operations. 
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Proof. We have by assumption that V (K) = S®Z where n(Z) = 0 and where the reduced 
echelon basis (t>i,..., v s ) of S is such that ..., tt(v s ) are non zero vectors in row 

echelon form. After reordering the columns as in the Proposition, it follows from the 
definition of the Hermite normal form that Tf{wj) = 7 which forces Wj —Vj € ZC |Z S . 
The conclusion then follows from Lemma 4.10. If the Wj’s have positive entries, the 
computation of all Fj’s reduces to a product of polynomials whose total degree sum is 
d y . This costs 0(d x d y ). □ 

The advantage is that we reduce our recombination problem to linear algebra (over 
Z): there are efficient algorithms to compute the Hermite normal form over Z with the 
same number of arithmetic operations as for the reduced echelon form over Q (26). 
The difficulty is that some of the 'Wj's might have negative entries, in which case the 
complexity of computing the Fj’s is more difficult to estimate. 

Example 4.14. Let us continue example 4.11. After reordering, we get that the Hermite 
normal form of the basis of V (K) D Z 4 is 

( Wl = ( 1 , 0 , 1 , 0 , 1 ), w 2 = ( 0 , 1 , 0 , 0 , 1 ),( 0 , 0 , 0 , 1 ,- 1 )), 

leading to the factorization 

(Fi,F 2 ) = ([F 1 F 3 F 5 ] d -+\ [F 2 F 5 ) d * +1 ). 

Here, the vectors wi and w 2 have positive entries, and the computation of the Fj’s is 
fast. Note that the vectors w\ and w 2 do not necessarily form a partition of (1,1,1,1,1) 
anymore. In particular, the analytic factor J -4 has disappeared from the recombination 
process due to the relation [F^[ dx+1 = [T^\ dx+1 . 

Example 4.15. Suppose that s = 5 and that 

s = ((1, 0 , 1 ,0,0), (0, 1 ,0, 1 ,1)) Z = ((0,0,1, 2 , -3). 

Then F is reasonnably ramified, the Hermite normal form of the basis of V (K) (~l Z 5 is 
((1,0,0,-2,3),(0,1,0,1,1),(0,0,1,2,-3)). 

The vector wi = (1,0, 0,-2,3) has now negative entries and the computation of the 
corresponding factor 



is a priori more expensive. Note that if we had reordered the indices such that Z = 
((0, 0, 2,1, —3), we would have obtained Wj and w 2 with positive entries. 

Remark 4.16. In general, we can show that if F is reasonnably ramified and if Z is 
generated by vectors with at most two non zero entries (meaning that all branches with 
high q -invariant intersect at most one other branch), then we necessarily have Wj £N S .A 
concrete example for which it is not the case is given by F = (y 6 — (y—x) 2 )(y—x) € Q[x, y]. 
In that case Z = ((0, 0, 0, 2, —1, —1)). 
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5. Conditions for constant residues. Equations of V(L). 

There are two main approaches that allow to determine when the residues of G/F do 
not depend on x. The first approach is related to the first De Rham cohomology group of 
the complementary set of the affine curve F = 0 in Aj|. It allows to test whether certain 
differential meromorphic forms are linear combinations of the logarithmic derivatives of 
the absolute factors of F by checking closedness. This is the approach followed by Gao 
(11). The second approach, that we will follow here, is based on a divisibility criterion by 
F and has been developped by Lecerf (20). Hence, this section is essentially a relecture 
of Section 1 in (20), except that we have now to take into account that the leading 
coefficient of F is not necessarily invertible in K[[x]]. 

From now on, we fix p. £ F s and we denote by G := G M the corresponding polynomial 
and by pi,. ■ ■, pa v the residues of G/F at the roots of F. We denote by p > 0 the 
characteristic of K and we adopt the convention K(x p ) = K for p = 0. 

5.1. Equations for E(K) 

We want to know when the residues pk £ K(x) lie in the subfield K. We recall that 
the usual K-derivation of K(x) uniquely extends to a K-derivation of K(a;). An obvious 
necessary condition for that p k £ K is that the derivatives of the p k s vanish. More 
precisely, we have the following lemma: 

Lemma 5.1. We have p' k = 0 if and only if pk £ K(x p ). If moreover K has characteristic 
zero or greater than 2 d x (d y — 1), then p' k = 0 if and only if pk £ K. 

Proof. See for instance the proof of Lemma 2.4 in (11). n 

Let us denote by K[x, the vector space of bivariate polynomials of degree < m 
in x and < n in y. We introduce the K-linear operator 


D : K[x,y\ dx ,d v -i 


G 


K[x, y]3d x -l,3dy~3 

(' G X Fy — GyFj^Fy ~ ^F X yFy ~ FyyF^G, 


with the standard notations F y , F xy , etc. for the partial derivatives. Let y k {x) be the 
root of F corresponding to the residue pk ■ By combining the formulas 


Pk(x) 


G{x,y k ) 
Fy(x,y k ) 


and y' k ( x) = 


F x {x,y k ) 
Fy{x,y k )' 


we are led to the equality 

_ D(G){x,y k ) 

X> Fy{ x > Vk) ’ 

so that 

Pk = 0 B(G)(x,y k )=0. 

In particular, since F is separable with respect to y, it follows that p' k = 0 for all 
k = 1,..., d y if and only if F divides D(G) in K(x)[y]. In order to reduce this division 
problem to a well estimated finite number of linear equations, we localize. 
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Let a € IK[x] be an irreducible polynomial which is coprime to the leading coefficient 
lc y (F) e K[cc] of F. We denote by 

A := IK[x] (a) 

the localization of K[x] at a. Hence, euclidean division by F in A[y] is well defined. Each 
Q € A[y] has a unique a-adic expansion 

+ OO 

Q(y) = ft £ K[a?, y], de gx ft < de g a - 

i =0 

For each pair of positive integers 0 < m < n, we introduce the truncated polynomial 

n— 1 

{oil ~ 

i=m 

We have the following lemma, generalizing Lemma 3 in (20): 


Lemma 5.2. Let F and D(G) as before and denote by D(G) = QF 
division o/D(G) by F in the ring A[y]. Let 


R the euclidean 


2 ( 4-1 

. deg a . 


+ 1 


and 


n := 


r3 d x — 1 
deg a 


1 . 


Then F divides D(G) in K(x)[j/] if and only if {Q}", = {l?}” = 0- 


Proof Since the leading coefficient of F is invertible in A, then F divises D(G) in K(a:)[y] 
if and only if it divises D(G) in A[y], Suppose that F divises D(G) in K(a;)[y]. Since both 
F and D(G) lie in IK[x, y] and F is primitive with respect to y, it follows from Gauss 
lemma that F divises D(G) in K[a;, y]. It follows that R = 0 and that Q e K[x,y] has 
degree deg^ D(G) — deg^, F < 2 d x — 1. In particular, the coefficients of a 1 in the a-adic 
expansion of Q are zeroes as soon as i deg a > 2 d x — 1, that is for all i > m. In particular 
{Q} m = 0. Conversely, suppose that {Q} m = {^?} = 0. Then the a-valuation of 

D(G) — {Q}" 1 A is greater or equal to n. It follows that either D(G) = {Q} m F or 

deg a ,(D(G) — |Q} m A) > ndega > 3 d x + a — 1. 

But we have that 

deg* {Q} m < (m — 1) deg a + deg a — 1 < 2 d x + deg a — 2 

so that deg a .(D(G) — {Q} m F) < 3 d x + deg a — 2. This forces equality D(G) = {Q} m F 
in A[y]. Since all members of the equality he in K(x)[y], the equality holds in the ring 
K(a;)[y]. □ 


Remark 5.3. If the leading coefficient of F does not vanish at 0, then we can take 
a(x ) = x and A = K[[a;]] in the previous lemma. More generally, if IK has cardinality 
greater than d y , we can reduce to that case up to replace F(x, y) by y dy F(x, a + l/y) for 
some a € IK such that F( 0, a) ^ 0. This Moebius transformation has a negligeable cost 
for our purpose. It only exchanges the points (0,oo) and (0, a) and does not modify the 
geometry of F along the fiber x = 0. 

If K has cardinality < d y , then both degenerate situations u(0) = 0 and F( 0, a) = 0 for 
all a € K might hold simultaneously. In such a case, we need to find a prime polynomial 
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a £ K[x] coprime to u. Note that we can find such an a with deg a £ G(log(deg u)) (see 
the proof of Proposition 6.1 in Section 4). 

Remark 5.4. The situation a(x) ^ x requires to perform euclidean division in A[y] with 
A ^ K[[x]]. This may break the sparse structure of D(G) inherent to the sparse structure 
of the analytic factors Jy £ K[[x]][j/]. 

According to Lemma 5.2, we introduce the following F-linear map: 

D a : F s —> K[x, y\ x K[x, y\ 



where a, m, n are defined as in Lemma 5.2 and where Q and s are defined by the euclidean 
division D(G M ) = QF + R in A [y]. 

Corollary 5.5. We have equality ker(D a ) = T4(]K(a: p )). If moreover K has characteristic 
zero or greater than 2d x (d y — 1), then ker(D a ) = V'(IK). 

Proof. Follows by combining Lemma 5.1 and Lemma 5.2. □ 

The image of D a is contained in the finite-dimensional vector space 

D a (F s ) c K[x,y]d x -i, 2 d y -3 x K[x,y] 3 d x -i,d y -i. 

Hence, for F = K, the computation of ker(D Q ) reduces to compute the kernel of a K-linear 
system of s unknowns and 0 {d x d y ) equations. 

5.2. The case of small characteristic 

When the characteristic p of IK is small, we need supplementary condition in order to 
know when pk £ K. In fact, we will get in such a case the stronger conditions pk £ F p 
thanks to the following F p -linear operator: 

N : K[x,y] dx}dy - 1 — > K[x, y p ] p d x ,d v -i 

G i—* G p — d^~ 1 (GF p ~ 1 ). 

The operator N is well defined since d y (N(G)) = 0. The vector space K[x a , 2 / fc ] mj „ has 
to be understood as the vector space of polynomials of bidegree (to, n) in the variables 
( x a ,y b ). The operator N was introduced by Niederreiter in the context of univariate 
factorization over finite fields (22), and then used for bivariate factorization in (20). 

Lemma 5.6. We have V(F P ) = {p £ F® | N(G M ) = 0} 

Proof. This follows from Theorem 2 in (22). □ 

Since deg a ,N(G p ) = pd x , the number of linear equations to be solved for computing 
ker(N(G^)) grows linearly with p. The idea developed in (20) is to combine N with the 
operator D in order to cut down this dependancy in p. 

Lemma 5.7. Suppose that p £ ker(D 0 ). Then N(G P ) £ K[x p ,y p ]. 
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Proof. See (20), Lemma 4. □ 

Hence, previous lemma ensures that the F p -linear map 

N a : ker(D a ) — >K[x p ,y p ]d x ,d v -i 
T'—> N(G /i ) 

is well-defined. In particular, if K = F p fe is a finite field, the computation of ker(D a ) 
reduces to compute the kernel of a K-linear system of dimker(D a ) < r unknowns and 
0(kd x d y ) equations over F p . 

6. Algorithms and complexity 

We combine now our previous results in order to give an algorithm for factorization 
and an algorithm for computing the number of rational factors, and we study their 
complexities. 

6.1. A rational factorization algorithm 

We obtain finally a deterministic algorithm for irreducible rational factorization of 
separable bivariate polynomials. The field F now stands for K if K has characteristic zero 
or greater or equal to d x {d y — 1) and F stands for the prime held F p of IK otherwise. 
Given a vector space V over F, we denote by rebF(H) the reduced echelon basis of V over 
F. If we say compute V, this means compute rebp(H). 

Algorithm : Critical Factorization 

Input: A bivariate polynomial F £ K[a;, y] separable with respect to y. 

Output: The irreducible rational factors of F. 

Step 0. Compute the content / £ K[x] of F with respect to y and do Ft— Fj f. Compute 
/i,..., ft the irreducible factors of / over K. 

Step 1. Compute the truncated analytic factors [Ji] dx+1 ,..., [■7 r s ] dx+1 of F. If s = 1 then 
return F is irreducible. Otherwise, build the polynomials \Fid y Fi\ dx+1 for all* = 1,..., s 
and initialize So 4— F s . 

Step 2. If ley (F)(0) ^ 0, let at—*. Otherwise, compute a £ F[*] of degree a £ 0(\ogd x ) 
such that a is irreducible and coprime to u. 

Step 3. Build the F-linear system associated to D a and compute So t— ker(D a ). If 
dimr So = 1, return F is irreducible. If K has characteristic zero or greater than d x {d v — 1), 
then go to Step 5. Else go to Step 4. 

Step 4. Build the F-linear system associated to N a and compute So t— So Oker(N a )). If 
dimp So = 1, return F is irreducible. If reb(S'o) does not form a partition of (1,..., 1) or 
if Z ^ 0 then go to Step 5. Otherwise, go to Step 7. 

Step 5. Compute q := \y/d\ where v := vala;(DisCy(F)) and d is the minimal y-degree 
of the [Fi} dx+1 . If q < d x , go to Step 7, otherwise go to Step 6. 
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Step 6. Compute the g-truncated analytic factors of F and compute So So D W q . 
Step 7. We have S = Sq. For each v 7 e reb(S), j = 1, ..., r, compute 


Fj := 


lc(F)n^T 


d x ~\~1 


and compute the primitive part Fj of Fj with respect to y. 
Step 8. Return /i,..., f t , F 1 ,..., F r . 


The proof of Theorem 1 follows from the following proposition. 

Proposition 6.1. The algorithm Critical Factorization is correct. It performs at most 
one univariate factorisation in K[.t] of degree d x and : 

- //IK has characteristic zero or greater than d x {d y — 1), then at most 

0{d x d 2 y + d v ma x(d x , q)s ul ~ 1 ) + C(max(<4, g)) 

arithmetic operations over K. 

- If IK = F p fe, then at most 

0(kd x d 2 + kd y ma x(d x , g)s“ _1 ) + (D(k)C(m&x(d x , q)) 
operations over F p . 

Proof. We have Sq = V(F P ) at the end of Step 4 by Lemma 5.1. Hence S = So if and 
only if Z = 0 by Proposition 4.1. If reb(So) does not form a partition of (1,..., 1) or 
Z yl 0, then S ^ So and we need to go to Step 5. Otherwise, the recombination problem 
is solved and we can go directly to Step 7. 

We have So = V(K) or So = P(F) at the beginning of Step 5. Since deg^, lc(J r i) < d x , 
we have that deg Fi = deg y [Fi] dx+1 . Hence q > N by the Remark 3.2. It follows from 
Theorem 3.4 and Corollary 3.5 that S = So at Step 7. For all j = 1,..., r, we have that 

Fj = since 

deg a .(k(F)/lc (Fj)) + deg X {F 3 ) < deg x F. 

Hence, the algorithm returns a correct answer. Let us study its complexity. 

Step 0. Computation of the content of F requires 0{d x d y ) arithmetic operations over 

K. 

Step 1. Computations of the ( d x + l)-truncated Ff s has complexity C(d x ) by defini¬ 
tion. Then we compute all \fiid y Fi\ dx+l with <D(sd x d y ) operations in K. 

Step 2. If IK has characteristic zero or greater than d x , we can take a = x — c for some 
c such that lc(F)(c) ^ 0. Otherwise, a basic approach (certainly not the most efficient) 
consists to remark that for n > log p (deg.j, lc(F)), the polynomial 

x pn — x 

gcd(x p " — x, lc(.F)) 

has positive degree and is coprime to u. Then, we take for a an irreducible factor of d n , 
which has necessary degree less or equal to n £ 0(logd x ). Thanks to fast gcd computa¬ 
tions, we compute a n within 0(max(deg x lc (F),p n )) C 0{d x ) operations. Then, we need 
to perform a univariate factorization of degree at most d x in order to find a. 
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Step 3. In order to build the linear system of the map D a , we need first to compute 
Di := D (G ei ) for e, varying over the canonical basis ofF s . This has complexity 0(sd x d y ). 
Then, we need to compute the a-adic expansion of the Df s and F. This costs 0(sd x d y ) 
thanks to (14), Theorem 9.15. We need then to perform s euclidean divisions in (A /a m )[y\ 
of polynomials of degree 0(d y ), and where m £ 0{d x /a). This costs 0(sd y ) operations 
in A /a m , hence 0(sd x d y ) operations over K. Then computing the reduced echelon basis 
of ker(D a ) requires 0{d x d y s u ~ 1 ) operations in K. 

Step 4. We compute here Sq with at most O^kdxdyS^ *) arithmetic operations over 
F p thanks to (20), Proposition 4. It has to be noticed that the construction of the linear 
system associated to N a might constitute a bottelneck to the algorithm. The F-vector 
space Z is determined by s unknowns and 0(d x d y ) equations over K, hence testing Z ^ 0 
requires at most 0(d x d y s u ~ 1 ) operations over K. Note that Z = 0 as soon as F is locally 
irreducible along x = 0 (Lemma 6.5). 

Step 5. Since v < 2 d x d y , we can compute Disc y (F) by considering Fas a polynomial 
in y with coefficient in the ring K[x]/ (x 2dxdy ). This requires 0(d x d y ) arithmetic opera¬ 
tions over K thanks to fast euclidean algorithm (14). Computing q then has a negligeable 
cost. 

Step 6. Computing the higher truncations of the Ffs requires C(N) operations over 
IK by definition. Then computing all [Fid y Fi]^ +1 requires 0{sNd y ) operations in K. 
Building the linear system that determines the equations of W N fl Sq requires at most 
0(Nd y ) operations in IK thanks to the sub-product tree technique. Then computing the 
reduced echelon basis of W N fl Sq requires at most 0((N — d x )d y t w_1 ) operations over 
K, where t = dimu U < s. 

Step 7. Let rij := deg y Fj. Then the computation of Fj requires 0(d x rij) operations in 
K, hence a total of 0{d x d y ) operations for all j. The cost of primitive parts computations 
amounts to the same number of operations. 

The proof follows by adding all these costs, and by remarking that one arithmetic 
operation over F p k requires 0(k) operations over F p . □ 

Remark 6.2. (About Step 5.) The cost of the discriminant computation might be very 
high in terms of bit complexity. Moreover, it might happen that N and q differ by a factor 
of the order of magnitude of d y . Hence it is much more preferable to approximate the qf s 
(hence N) during the computation of the FC s, for instance by using the relations with 
the characteristic Puiseux exponents of the branches Fj (see for instance (23)). Hence, 
the cost of Step 5 is included in C(N). Note that if the reduced echelon basis of V(K) 
does not form a partition of (1,..., 1) or if Z ^ 0, then necessarily N > d x . Finally, note 
that if IK = Q (or more generally a number field), there are strategies to compute the 
valuation of the discriminant by working modulo a well chosen prime p (23). 

6.2. Algorithms for the number of irreducible rational factors 

In the first algorithm, we suppose that IK has positive characteristic, or that K is a 
decomposition field of F. The field F stands for K if K has characteristic zero or for the 
prime field F p of K otherwise. 

Algorithm : Number of Factors 

Input: A bivariate polynomial F £ K[x, y] separable with respect to y. 
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Output: The number of irreducible rational factors. 

Step 0. Compute the content / £ K[x] of F with respect to y and do F <— F/f. Compute 
t the number of irreducible factors of / over K. 

Step 1. Compute the d x + 1-truncated analytic factors of F. 

Step 2. Compute r £- dim V(F). If r = 1 then return t + 1. 

Step 3. Compute r £- r — dim Z. Return t + r. 


Proposition 6.3. The algorithm Number of Factors is correct. It performs at most one 
univariate factorisation in IK [a;] of degree d x and : 

- //K has characteristic zero, then 0(d x d y s u ~ 1 ) +C(d x ) operations over K. 

- If K = F p fc, then 0(kd x d y s ul,_1 ) + 0{k)C(d x ) operations over F p . 

Proof. The correctness of the algorithm follows from Proposition 4.1. In positive charac¬ 
teristic, we can compute P(F) as in the previous algorithm and the complexity analysis 
follows from the proof of previous proposition 6.1. In characteristic zero, we have that 
V(F) = V (K) by Lemma 4.3 so that we compute V (F) as in Step 3 of algorithm Critical 
Factorization. The F-vector space Z is determined by s unknowns and 0(d x d y ) equations 
over K, hence testing Z ^ 0 requires at most 0(d x d y s ul ~ 1 ) operations over K. □ 

Finally, we have the following algorithm for an irreducibility test. Here, no hypothesis 
are made on the field IK. 

Algorithm : Irreducibility Test 

Input: A bivariate polynomial F £ K[x, y] separable with respect to y. 

Output: True if F is irreducible over IK, False otherwise. 

Step 0. If F is not primitive with respect to y , then return False. Otherwise, replace F 
by its primitive part. 

Step 1. Compute the n-truncated analytic factors of F with n = d x + 1 if Char(K) > 0 
and n = 2 d x otherwise. 

Step 2. Compute r <— dimH(F) if Char(K) > 0 and r <— dimV r (K) D W 2dx otherwise. 
If r = 1 return True. 

Step 3. Compute r •<— r — dim Z if Char(K) > 0 and r •<— r — dim Z n W 2dx otherwise. 
Step 4. Return True if r = 1, False otherwise. 


Proposition 6.4. The algorithm Irreducibility Test is correct. It performs at most one 
univariate factorisation in K[x] of degree d x and : 

-If K has characteristic zero, then 0(d x d y s u ~ 1 ) +C(2d x ) operations over¥L. 

- If K = F p fc, then 0(kd x d y s u> ^ 1 ) + 0(k)C{d x ) operations over F p . 
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Proof. The correctness of the algorithm follows from Proposition 4.1 in positive char¬ 
acteristic and from Proposition 4.7 otherwise. The complexity analysis follows from the 
proof of the two previous algorithms. □ 

The proof of Theorem 3 follows from the two previous propositions 6.3 and 6.4. 

6.3. Locally irreducible polynomials 

We recall from the introduction that we say that F is locally irreducible along the 
line x = 0 (resp. absolutely locally irreducible) if the germs of curves ( C,P ) C (Pj|,P) 
defined by F are irreducible over K (resp. over K) at each rational place P of the line 
x = 0, including the place at infinity. 

Lemma 6.5. Suppose that K has characteristic greater or equal to d y . If F is locally 
irreducible along the line x = 0, then Z = 0. 

Proof. Let p G Z. we have in particular 

S 

Y / p i F i d y F i (0,y) = 0. (11) 

i =1 

Let us consider first an analytic factor Fj corresponding to affine place of C along x = 0. 
Since F is locally irreducible along the line x = 0. It follows from Hensel’s lemma (14) that 
Fi(0,y) is a power of a prime polynomial that is coprime to Fi(0,y). Hence relation (11) 
combined with Gauss Lemma imposes that Fi(0,y) divises pid y Fi(0,y), hence Pi = 0 
thanks to the assumption on the characteristic of K. Suppose now that F has an analytic 
factor, say T\ that vanishes at (0,oo). By the local irreducibility assumption, F\ is the 
unique such factor. Hence p € Z becomes equivalent to that 

P 2 = ■ ■ ■ = p s = 0 and p 1 [f 1 d y (Fi)] dx+1 = 0, 

thanks to what we proved for the affine places. The leading coefficient of T\ d y {F\ ) is 
equal to dilCy(F). It has degree < d x , and rfi ^ 0 by assumption on K. Hence, last 
equation implies p± = 0. □ 

Lemma 6.6. If Card(K) > d y and F is locally irreducible along the line x = 0, then we 
compute the (d x + 1 )-truncated analytic factors of F with one univariate factorization of 
degree at most d y plus 0(d x d y ) arithmetic operations over K. 

Proof. Since Card(K) > dy, there exists yo £ K such that F(0,y o ) ^ 0. To find such an 
element yo has a negligeable cost once the univariate factorization of F( 0, y) is given. Then 
the Moebius transformation F «— y dy F(x,a + 1 /y) reduce to the case where u = lc V (F) 
is a unit modulo x. Hence we are in position where 

F(0, y ) = u(0)f[pr, Fi(p,y) = pr 

i= 1 

for some distinct prime polynomials Pi € K[y]. We can lift this factorization modulo x dx+1 
with 0{d x d y ) arithmetic operations over K thanks to the multifactor Hensel lifting (14), 
Theorem 15.18. Then, we perform the inverse Moebius transformation in order to get the 
analytic factors of the original polynomial. Both Moebius transformations can be done 
in softly optimal time 0(d x d y ) with interpolation/evaluation. □ 
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Proof of Theorem 2. It follows from Lemma 6.5 and Proposition 4.1 that S = V (IK) under 
the assumption of Theorem 2. If K = F p k is a finite field, we can compute a basis of V (K) 
from the (dr + l^truncated factors within 0{kd x d v s u> ~ 2 ) operations over F (see the proof 
of Proposition 6.1). Since one operation in IK amounts to 0(K) operations over F, the 
proof follows from Lemma 6.6. If IK is a decomposition held for F, we have V (IK) = V (K) 
thanks to Lemma 4.3. Since IK has characteristic zero or greater than d x (2d y — 1), we 
compute the reduced echelon basis of V(K) within 0{d x d v s ul ~ 1 ) operations over K. □ 

7. Absolute factorization 

In order to generalize our results to the absolute case, we follow the strategy developed 
by Cheze-Lecerf in the regular case (9). 

Absolute factorization. We denote by E\,... ,E? the irreducible factors of F in IK [a;, y]. 
We represent this factorization by a family of pairs of polynomials 

{(Pi,<7i),...,(P t ,g t )} 

where qk G IK[z] is monic, deg z Pk < deg qk, Pk(x,y,<p) G IK[x, y] has constant bidegree 
when <f> runs over the roots of qk, and for each factor Ei there exists a unique pair (fc, </>) 
such that qk{4>) = 0 and 

Ei(x,y) = Pk{x, y, (f>). 

Such a representation is not unique, but is not redundant. The qk s are irreducible if and 
only if the products n<j t (0)=o -Pfe(^) are the irreducible factors of F in K[x, y]. 

Absolute analytic factorization. We denote by £i,... ,£§ the irreducible analytic factors 
of F in K[[x]][y]. As before, we suppose that the £fs are given by a collection of pairs of 
polynomials 

{CPi,pi),...,{Ve,pe}} 

where pk G IK[z] is monic, deg z Vk < degpfc, Vk{4>) £ K[[a;]][y] has constant degree in y 
when <j) runs over the roots of pk, and for each £i there is a unique pair (k, <f) such that 
Pk(<t> ) = 0 and 

£i{x,y) = Vk{x, y, </>). 

In particular, the pk s are separable. The pk s are irreducible if and only if £ = s, if and 
only if 

(deg y (Vk),deg z (pk)) = (e fe ,/ fc ), 

where e*, and fk stand for the ramification index and residual degree at the rational 
places of F = 0 over x = 0. We do not necessarily assume this here, the only important 
point from a complexity point of view being that we necessarily have 

i S 

8 = = H/i- 

k— 1 i= 1 

In particular, the more the curve F = 0 is ramified over x = 0, the smaller the number 
of unknowns is. This is a great difference with the regular case, for which equality s = d y 
always holds. We call the n-truncated absolute analytic factorization the data of the pairs 
([ V k ] n+ \Pk)- 
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Solving recombinations with K -linear algebra. In analogy to the rational case, we denote 

by 

S = (tii,.. .,u f >K C K s 

the K-vector space generated by the recombination vectors Vi,... ,v 3 solution to 

S 

e j ='«j n y ’ 

i=i 

with Uj £ K[x], Uj( 0) = 1. For /i G I s , we denote by 

s 

Qfi ■= Y, A ijEjdyEj G K[[x]][y]. 

i=l 

We introduce the K-vector spaces 

V := {/X G K* I [g^ +1 G , ErdyE?)}. 

Hence p £ V if and only if the residues of [Q^\ d:c+1 /F lie in K by Lemma 3.8. We 
introduce also 

Z~{p£ K s ‘| [Gfj\ dx+1 = 0}, 

and 

W := (m G k s I [g^l +1 = 0} 

where N is the separability order of F. Lemma 3.8 combined with Theorem 3.4 and 
Proposition 4.1 give the relations 

S = VnW and V = S ® Z. 

First equality leads to an algorithm for solving recombinations. The second equality leads 
to an algorithm for computing the number of irreducible absolute factors. The idea now 
is to use the Vandermonde matrices attached the polynomials p k s in order to compute 
a basis of the involved vector spaces with linear algebra over K. 

The Vandermonde isomorphism. We define the partition {1,..., s} = 7i U • • • U I( by 
requiring that 

£i(x,y)= JJ Vk{x, y, <f>), k = 

i£lk Pk(4>) = o 

These products lie in K[[x]][y]. They coincide with the irreducible analytic factors of F if 
and only if the p k s are irreducible. Let nk := degpfc. If p G K s , we denote by G IK"* 
the vector whose entries are the entries of p whose index lie in Ik- We introduce the 
K-linear map 

A = (Ai ,..., At) : K* —> ®k K 

p = {p (1) ,... ,p W ) i— >v:= . ,u (e) ) 

where Ak stands for the transposed of the Vandermonde matrix of the roots (<pki , • • •, 4>kn k ) 
of pk- In other words, the vector 

P k) := A k pW G K Jlfe <g> K 

is defined by 

nk 

■f ■■= X>f C 1 - 
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Since the polynomials pk are separable, each map Ak : K" fe —► is an isomorphism, 

hence so is the map A. 

Recombinations over K. For a given Q £ K[[a;]][y, z], we denote by coeff(<7, z £ K[[a;]][y] 
the coefficient of zi in Q. We have by construction that Vk divises F in the ring I€[.s]/(Pk)[[®]][y] 
and we denote by Vk the unique polynomial such that F = VkVk £ K[[x]][2/, z], with 
deg z Vk < deg pk- Given v £ K s , we denote by 

t jik 

= EE coeff (VkdyVk ,z J 1 ). 

fe=i i=i 

We introduce the K-vector spaces 

:= [u £ K* | [H v ] d ’ +1 £ (EidyEi ,..., E f d y E f ) K } 

Z K := |^£K s " | [H„] d * +1 =0}. 

and 

W k ~{v£K b | [W„]£ +1 = 0}, 

We have the following 

Proposition 7.1. The isomorphism A : p K > v induces isomorphisms 
V = Vk K, W = Wk 8>k K and 2 = Zk K. 

Proof. By construction, we have that 

t rife i n k n k 

0, = £E P^VkdyVkiM = EEE coeff (VkdyVk, 

k—1 i—1 k—1 i= 1 j=l 

^ rife rife 

= EE(E Ki X ) coeS(Vkd y Vk,z J *) 

/C=l J = 1 2—1 

L n k 

= EE ^ (fe) coeSiAdyVk,^- 1 ) =H„. 

k= 1 J=1 

The claimed isomorphisms then follow from the definitions of the involved vector spaces. □ 

Proof of Theorem f. We have shown how to compute a basis of all involved vector spaces 
with linear algebra over K. Unfortunately, we don’t have the reduced echelon basis trick 
when working with the unknowns v instead of p. To solve this problem, we rather use an 
absolute partial fraction decomposition algorithm along a regular fiber, following Section 
4 in (9). We obtain the following algorithm. 

Algorithm : Absolute Factorization 

Input: A field K with cardinality at least d x (2d y — 1) and a bivariate polynomial F £ 

K[x, y] separable with respect to y. 

Output: A family {(Pi, qi ),..., (Pj, qt)} that represents the absolute factorization of F. 
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Step 1. Compute a basis u-[,..., v f of tic D Wr. 

Step 2. Find a regular fiber x = a for some a € K. 

Step 3. Compute hi := TL V1 (a, y ),..., hf := TL^{a , y). 

Step 4. Call Algorithm 7 in (9) in order to find (ci,...,Cf) € K r that separate the 
residues of the hf s. 

Step 5. Let h = Cj/ij. Call the Lazard-Rioboo-Trager algorithm (Algorithm 14 in (9)) 
in order to compute the absolute partial fraction decomposition of h/F(a,y). 

Step 6. Call Algorithm 6 in (9) of absolute multi-factor Hensel lifting in order to lift the 
decomposition of Step 6 to {(Pi, <71), ... ,(P t ,q t )}. 

Proposition 7.2. Let m = ma x(d x + 1 ,N) and p = Char(K). The algorithm Absolute 
Factorization is correct. It performs at most 

ma x(d x + 1, N)d y + rd x d y ) + C(max(d x + 1, N)) 

operations mKj/p^O or p > dx(2d y — 1) and at most 

0(ks UJ ^ 1 max(da; + 1, N)d y + kfd x dy) + 0(k)C(ma,x(d x + 1, N)) 

operations over F p if K = F p fc. 

Proof. Given the 7Vs and the pk s, we can compute a basis of the K-vector spaces Vr 
and W K with the the same cost as in the rational case, with the number of unknowns s 
being replaced by s. Given {y\,... ,Vf) a basis of Vr fl Wr, we have by construction that 

(’Hvi , • • ■ , Pvrjt. = {EldyEl, . . . EfdyEf) R, 

with TL Vi € K[x, y]. By assumption on the cardinality of the field, we know that there 
exists a regular fiber x = a over which F(a,y ) is separable of degree d y . We can find 
such a fiber by computing Disc F(i, y) for i = 1,..., d x (2d y — 1) + 1 until we reach a 
non vanishing discriminant. This costs at most 0(d x d y ) operations over K. Then we 
refer to (9), Paragraph 4 for the remaining steps of the algorithm. Step 4 costs 0(fd x d y ) 
operations in K, step 5 costs 0(d y ) operations in K and step 6 costs 0(d x d y ) operations 
in K. □ 

Theorem 4 follows immediately from the previous proposition. 

Proof of Theorem 5. We have by Proposition 7.1 that the number of absolutely irre¬ 
ducible factors of F is equal to 

f = dimjj S = dim K Vk ~ dim K Z K 

Given the d x + 1-truncated analyic factorization of P, we can compute dimR Vr and 
diniR Zr with the same costs as in the rational case, with the number of unknowns s 
being replaced by s. The proof of Theorem 5 follows. □ 
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8. Non degenerate polynomials 


We introduce the notion of P-adic Newton polytopes. These combinatorial objects 
give a lot of interesting informations for both rational and analytic factorization. In 
particular, we show here that they permit to detect a large class of polynomials whose 
separability order is small. 

Let us fix P G K[y] a non constant polynomial. Any polynomial T G K[[:r]][y] can be 
uniquely expanded as 

P(x, y) = ^2 fijX'P’ G K[[x]][y], 
with fij G K[y], deg fij < degP. Let 

Suppp(P) := {(i,j) G N 2 , f i:j ± 0} 

stands for the P-support of T. The P-adic Newton polytope of P, or P-polytope for 
short, is the convex hull of the positive cone generated by the support of T, that is 

N p : jr ■= Conv ^(Suppp(P) + (R + ) 2 ^. 

When P = y we recover the usual notion of Newton polytope of a bivariate power series 
(17), and we might simply say Newton polytope for the y-polytope. Take care that the 
terminology of Newton polytope refers sometimes in the litterature for the (compact) 
convex hull of the support of a bivariate polynomial. We have the following lemma. 

Lemma 8.1. Let P G K[y] be separable and irreducible and let a G K be a root of P. 
Then, the P-polytope of T coincides with the Newton polytope of F(x,y + a). 

Proof. Since P is irreducible and separable, the highest power of P that divises a given 
/ G K[y] coincides with the highest power of y — a that divises / in K[y], which coincides 
with the highest power of y that divises /(y + a). The Lemma follows. □ 

We call the P-edges of T the compact edges of its P-polytope. Let A be a P-edge and 
let oa and b\ stand respectively for the distance from A to the y-axis and a;-axis. We 
define the P-edge polynomial of F associated to a A as 

fp,A := x~ aA y~ bA ^2 fijx'y 3 e Kp^y]- 
A 

where fij G Kp := K[y]/(P) stands for the reduction modulo P. By construction, the 
polynomial fp_A is quasi-homogeneous and monic with respect to x and y. We say that 
a series is P-convenient if it is not divisible by P or x. 

Definition 8.2. We say that T G K[[x]][y] is non P-degenerate if it is P-convenient and 
if both P and all the P-edges polynomials of T are separable with respect to y. We say 
that T is non degenerate at infinity if y dy F{l/y) is non y-degenerate. We say that T is 
non degenerate if it is non P-degenerate for all irreducible factors P of P(0,y) and if it 
is non degenerate at infinity. 

Remark 8.3. By quasi-homogeneity, we can let x = 1 for checking separability of the 
P-edge polynomials. 
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Remark 8.4. Usually, the notion of non degenerate polynomials in K[[x]][y] only allows 
ramification at the places y = 0 and y = oo, while we consider here all places of Pjjj. A 
notable exception is (24) where the authors use collection of P-adic polytopes in order 
to improve the usual Bernstein-Koushnirenko bound for the number of solutions of a 
polynomial system with isolated roots. 

Remark 8.5. In zero characteristic, non y-degeneracy is equivalent to the most common 
definition of non degeneracy introduced by Kouchnirenko (17). In positive characteristic, 
there are several notion of non degeneracy. Kouchnirenko non degeneracy is equivalent 
to that the edge polynomials are separable with respect to x and y. This is the one that 
allows to generalize the Milnor formula to positive characteristic. Our notion is weaker 
(for instance y 3 — x 2 in characteristic 2). Weak non degeneracy introduced in (6), Section 
3 is equivalent to that the edge polynomials are squarefree. This is the one that allows 
to compute the number of local factors. Our notion is stronger (for instance y 3 — x 2 in 
characteristic 3). 

Lemma 8.6. Let P £ K[y] be separable and irreducible. Then F is non P-degenerate if 
and only if F(y + a) is non y-degenerate at any roots a of P. 

Proof. Let a be a root of P and let us write F(y + a) = Y^ c ij xl y :i f° r some dj £ K. 
A straightforward computation shows that the coefficients in the two expressions are 
related by 

Cij = P'(ay fij (a). 

By Lemma 8.1, the y-edges of F(y + a) are one-to-one with the P-edges of F. Let A 
be such an edge. The corresponding y-edge polynomial f yi A of F(y + a) and P-edge 
polyonomial f p a of J- are related by the formula 

fy,A ( X ’ pqK) = x ~ aA y~ bA H fij(a)x l y J = ev Q (ip, a) e K(a)[x,y]. 

[a> (i,j)e A 

where ev Q is induced be the isomorphism Kp ~ K(a) determined by a. Since the dis¬ 
criminant of monic polynomials commutes with specialization, we deduce that f y ,A is 
separable with respect to y if and only if fp ; a is. □ 

Proposition 8.7. Suppose that F £ K[x, y] is non degenerate, then S = V(K). 

Proof. Let as usual F = uF^ ■ ■ ■ F s be the irreducible factorization of F in K[[x]][y]. 
By Corollary 3.5, it’s enough to prove that qi < d x for all i. By point (2) in Lemma 
3.1, the y-invariant of the irreducible factors of T r in K[[x]][y] are all equal to yp Hence, 
there is no less to suppose that K = K. In such a case, Hensel lemma implies that we 
necessarily have P,;(0,y) = (y — oy) di for some on £ IK. Since qi is invariant under the 
Moebius transformations F ■£- P(x,y + af) (or F £- y dy F(x, 1/y) for ay = oo), we are 
reduced to estimate y; when P)(0, 0) = 0. Clearly, y, only depends on those factors Fj for 
which Pj(0,0) = 0 so that we can suppose that P £ K[x, y] is a product of Weierstrass 
polynomials which by Lemma 8.6 is non y-degenerate. By the multiplicative property of 
the resultant, we have 

dtqi —val„ (Disc, y ( F ,)) T ^ val 7: H.esy (F /, Fj ). 
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For given two Weirestrass polynomials J-,Q& K[[x]][?/], we have the formula 

valz ReSjX-F, Q) = (P, Q ) 0 


where (F,G)o stands for the intersection multiplicity 


(P ,£)0 := dim K 


i?,S) 


of T and Q at (0,0), see for instance (27) p.28 (the proof adapts to the positive charac¬ 
teristic case). Let us denote by 


A* := (R + ) 2 \ Njr. 


the complementary set in (K+) 2 of the Newton polytope of p. Note that Aj is compact 
since p is convenient by assumption. By Bernstein-Khovanskii-Kouchnirenko theorem, 
we have the formula 

(Pj,P,)0 > [Aj, A j], 

with equality if the product TiTj is non y-degenerate (the converse holds in characteristic 
zero). See for instance Corollary 5.6 in (8) in the case IK = C or (17) for any algebraically 
closed field. Here, 


[A i; Aj\ := Vol(Aj + Aj) - Vol(Aj) - Vol(A j) 
stands for the mixed volume of polytopes. In the same way, we have that 

(Pi, d y Ti )0 > [Aj, Aj] - d v , 

with equality if Pj is non degenerate (see Theorem 5.6 in (8), the proof adapts to the 
positive characteristic case since fv is assumed to be separable with respect to y). Since 
all TVs are irreducible, it follows that Aj is a triangle with vertices (0,0), (aj,0), (0,dj) 
where aj = val x (Pi (ir, 0)) (with aj = 0 and Aj being a segment if Ti does not depend on 
x). In such a case, we get that 

[Aj, Aj] = mm{diCLj,djai}, 

see (8), Section 5. Hence it follows that 

S 

qidi < ^2 didj - d y = di val x F(x, 0 ) — d y < did x - d y . 

3 =1 

The inequality <?j < d x follows. □ 


Let F G K[x, y] and suppose given the irreducible factorization 

t 

F(o, y ) = i[pr 

i=1 

For each i, we denote by Sj the total number of irreducible rational factors of the Pj-edges 
polynomials of F and by li the lattice length of the Pj-boundary. Note the inequalities 

Sj < ii < Hi. 


In the same way, denote by Soo the total number of rational irreducible factors of the 
edge polynomials of y dy F(x, 1/y) and by i oo the lattice length of its Newton boundary. 
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Lemma 8.8. Suppose given F £ K[x,y] primitive with respect to y and x. Then, the 
respective numbers s and s of irreducible analytic factors of F over IK and IK satisfy 

t t 

s < SF := 3j + Sqq and s < s F := ^ tj deg Pj + 4o, 

2=1 2=1 

both inequality being equalities if F is non degenerate along the fiber x = 0. 

Proof. Suppose first F monic w.r.t y. By Hensel Lemma, the number of irreducible 
factors of F in K[[x]][y] is the sum of the numbers of irreducible factors in the local 
rings K[[x]][y](p ; ) when Pj runs over the irreducible factors of P(0, y). Then the assertion 
for s follows for instance from Chapter 6 in (7). For the absolute case, we consider the 
decomposition 

S[N][2/](Pi) = 0 K[[x,y-a]]. 

P(a)=0 

By Lemma 4.10 in (6), the number of absolute factors of F in JK[[x]][y]( y _ a ) is bounded by 
the lattice length of the Newton boundary of F(x,y + a), which by Lemma 8.1 coincides 
with £i. Moreover, there is equality if P(x, y + a) is non y-degenerate, which is the case 
when F is non degenerate by Lemma 8.6. Summing up over all the roots of Pi and over 
all i, we get the result. If F is not monic, we conclude in the same way by taking also 
into account the place at infinity, n 

Remark 8.9. Equalities s = sf and s = sp hold in fact with the weaker hypothesis 
that the edge polynomials are squarefree. 

Example 8.10. Suppose that 

F(x,y) = y(y 2 - 2) 3 - x 2 (y 2 - 2) + x 5 £ Q[x,y]. 

Then F(0,y) = y(y 2 — 2) 2 has two irreducible coprime factors Pi = y and P 2 = y 2 — 2. 
There is a unique Pi-edge polynomial f'A.p, = —8 y + 2x 2 , which is obviously separable 
and irreducible. Hence s 1 = £1 = 1. There are two P 2 -edge polynomials 

fi := </>y 2 — x 2 and f 2 := — y + x 3 

where cj> £ Qp 2 stands for the residue class of y. Both polynomials are separable. Since <j> 
is not a square, fi is irreducible over Qp 2 , but has two factors over Q. Obviously, f 2 has 
exactly 1 factor over any field. Hence S 2 = 1 + 1 = 2 while £2 = 2 + 1 = 3. Since there 
are no points at infinity, we finally get that F has exactly s = si + S 2 = 3 irreducible 
factors in Q[[x]][y] and s = £ 1 + 2£ 2 = 7 irreducible factors in Q[[x]][y]. 

Example 8.11. Suppose that 

F(x, y) = y 6 (y 2 + l) 15 - x 10 (l + y 21 ) £ K[x, y], 
where IK is any field of characteristic 2. Then 

F(0,y) = y 6 (y 2 + l) 15 

has only two distinct irreducible factors P± = y and P 2 = y 2 + 1. There is a unique Pl¬ 
edge polynomial hyp, = y 6 — x 10 . It is separable with respect to y if and only if p ^ 2, 3, 
and we have si = £\ = 2 if p ^ 2. There is a unique P 2 -edge polynomial 

fA ,p 2 = </>V 5 - x w (l + </> 21 ) = -y 15 + (1 + 4>)x w , 
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where <p stands for the residue class of y in Kp 2 = K[y]/(y 2 + 1). This polynomial is 
separable if and only if p ^ 2,3, 5 and it is squarefree if and only if p ^ 2. For p ^ 2,3, 5, 
it has S 2 = 2 irreducible factors over Kp 2 while the lattice length of A is l 2 = 5. Hence, 
if p ^ 2, 3, 5, then F is non degenerate and has s = si + S 2 = 4 factors in K[[a;]][j/] and 
s = £ 1 + 2£ 2 = 12 factors in K[[x]][y]. 

Corollary 8.12. There is a deterministic algorithm that, given F £ K[x, y], returns 
False if F is degenerate and returns the irreducible factors of F over K otherwise. The 
cost of this algorithm is one univariate factorization of degree at most d y and 

0(d x d y a^- 1 )+C(d x ) 

operations over K i/p> d x (2d y — 1) or 

OikdxdyS"- 1 ) + 0(k)C{d X ) 
operations over F p if K = . 

Proof. The algorithm is as follows. We first compute the factorization F(0, y) = n!=i -F* n< ■ 
For each i, we test the separability of Pi and of all the Pledges polynomials of F. This 
step costs at most 0(d x d y ) operations over K. If F is non degenerate along the fiber 
x = 0, then s = sf by Lemma 8.8 and the separability order satisfies N < d x by 
Proposition 8.7. Hence we can take q = d x in Theorem 1. Corollary 8.12 follows. □ 

Remark 8.13. Analytic factorization of non degenerate polynomials may be reduced to 
Newton iteration in the rings Kp[[x]][y] after some translations and monomial change of 
coordinates. A first estimate leads in that case to C(d x ) C 0(spd x d 2 ), but we believe we 
can do better. One of the main obstruction is the difficulty to use a dichotomic multifactor 
Hensel lifting as in the regular case. 

Remark 8.14. A cheap pretreatment of F is to look at the fibers x = 0 and x = 00 (or 
y = 0 and y = 00 by reversing the roles played by x and y) in order to check if there 
is a fiber over which F is non degenerate and with the smallest sp or sp as possible. 
If we take for instance F(x,y) = y 6 (y — l) 15 — x 10 + x 9 y 21 , then the fiber x = 0 leads 
to Sp = 4 and sp = 7 while the fiber x = 00 would lead to s = s = 1 from which we 
immediately deduce that F is absolutely irreducible, whatever the field is. This kind of 
strategies based on the relations beewteen the (global) Newton polytope and bivariate 
factorization have already been considered in the litterature, see for instance (12) or (29) 
and the references therein. 


9. Conclusion 

When compared to the regular case, a great advantage of working along a critical 
fiber for factorization is that one has in general less analytic factors to recombine (always 
strictly less in the absolute case K = IK). Unfortunately, this might require in general a 
higher truncated precision, and our main Theorem 3 seems to suggest a bad worst case 
complexity. However, the important classes of non degenerate polynomials and locally 
irreducible polynomials illustrate that we sometimes gain in complexity when compared 
to the regular case. Might this hold in all generality ? We discuss briefly the two main 
obstructions for this. 
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9.1. Fast analytic factorization ? 

The strength of our approach deeply relies on the complexity of ./V-truncated analytic 
factorization. This is a crucial problem in singularity theory. One approach in charac¬ 
teristic zero or > d y is to compute the rational Puiseux series. There are well known 
algorithms for this, with actual complexity 0(d 5 ) in terms of the total degree d of F 
in the case of finite fields (23). This complexity is a bit too large for our purpose, but 
the authors told us that they recently developped a 0(d 4 ) complexity algorithm. This 
is the same complexity as for absolute factorization. In order to fit also in the rational 
factorization complexity class, we would need 0(d u+1 ) for analytic factorization. This is 
an open problem. 

Note that there exists algorithms for testing local analytic irreducibility of a germ 
of curve that do not require the use of Puiseux series. The main ingredient is that of 
approximate roots and makes essentialy use of resultants computations (1), (19), (10). 
Up to our knowledge, no complexity analysis have been done yet. We don’t know if 
approximate roots could lead to fast analytic factorization that would avoid Puiseux 
series computations. 

In characteristic p < d v , the concept of Puiseux series does not make sense and analytic 
factorization is a much more delicate problem. See for instance (15) for the generalization 
of Puiseux series in small positive characteristic. 

9.2. Fast recombinations ? 

Even if we get a fast algorithm for computing the iV-truncated analytic factors, our 
brute force analysis of the recombination problem in this paper led to a complexity 
0(Nd y r ul ~ 2 ) C 0(d ul+2 ). However, the polynomials [Fid y F,] V+1 we want to recombine 
have a very particular structure. Namely, they vanish with very high order at some points 
of the curve F = 0. An approach could be to write these polynomials in a basis constitued 
of adjoint polynomials. The number of required equations for solving recombinations 
(divisibility test by F or closedness of differential forms) would then decreaze. This fact 
is illutrated in some of our previous works, where we studied the relations beetween 
resolution of singularities and factorization (30) or toric genometry and factorization 
(29). Namely we developed factorization algorithms whose linear algebra steps belongs 
to 0{jp a s u J ~ 1 ), with p a < d 2 being the arithmetic genus of the strict transform of the 
curve F = 0 after some sequences of blowing-ups. 
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